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Spectroscopy of vortex bound states can provide valuable information on the structure of the 
superconducting order parameter. Quasiparticle wavefunctions are expected to leak out in the 
directions of gap minima or nodes, if they exist, and scanning tunneling spectroscopy (STS) on 
these low-energy states should probe the momentum dependence of the gap. Anisotropy can also 
arise from band structure effects, however. We perform a quasiclassical calculation of the density 
of states of a single vortex in an anisotropic superconductor, and show that if the gap itself is 
not highly anisotropic, the Fermi surface anisotropy dominates, preventing direct observation of 
superconducting gap features. This serves as a cautionary message for the analysis of STS data on 
the vortex state on Fe-based superconductors, in particular LiFeAs, which we treat explicitly. 

PACS numbers: Valid PACS appear here 



Introduction. Four years after the discovery of the iron- 
based [11 [2] high temperature superconductors, the struc- 
ture and symmetry of the gap function are still being de- 
bated. There is considerable experimental evidence that 
there is no universal gap shape [SHS] , perhaps in part due 
to the electronic structure that combines small electron 
and hole pockets, leading to an "intrinsic sensitivity" |6] 
to details. It is likely [5] that in most cases the gap has 
Aig symmetry, which, however, allows a continuous de- 
formation from a full gap to that with nodes on the Fermi 
surface (FS) sheets. Bulk experimental probes of the gap 
structure include specific heat and thermal conductivity 
oscillations in an external magnetic field [7, 8 , performed 
on the Fe(Te,Se) system [9 and P-doped 122 family [10] 
respectively. In both systems the oscillation pattern was 
found to be consistent [8, 10-12 with an anisotropic gap 
with minima along the F — X axis (in the unfolded Bril- 
louin zone), as predicted by spin fluctuation theories (see, 
e.g. Ref. 5). 

Order parameter structure is also reflected in the lo- 
cal properties of inhomogeneous superconducting states. 
Inhomogeneities may arise due to impurities, and the re- 
sulting quasi-bound states in nodal superconductors have 
tails that "leak out" in the nodal directions [T3J, provid- 
ing a signature of the amplitude modulation of the gap. 
The interpretation of these impurity states is complex: 
disorder potentials can be of the order of electron volts, 
and hence relatively high energy processes control the 
formation of such states, as well as their contribution to 
scanning tunneling spectroscopy (STS) images [14 . 

Under an applied magnetic field, inhomogeneous su- 
perconductivity arises due to modulation of the order 
parameter in a vortex lattice, and bound states localized 
around the vortex cores appear. In this case, relevant 
energy scales are of the order of the gap or lower and 
the bound states properties are determined by the shape 
of the gap and the band features near the Fermi sur- 



face. The decay length of the core states is of order of 
^0 = vf/ttA^ where vf the Fermi velocity and A is the 
gap amplitude. Consequently, variation of the gap with 
direction k at the FS, A(k) ^ const, directly influences 
the shape of the core states in real space, leading to the 
"tails" extending along nodes or minima. Since the de- 
cay of these states is exponential in distance p from the 
center of the vortex (except along true nodes where it 
follows power laws), these tails are very clearly seen in 
local measurements, and can be used to probe the gap 
shape [15 . Difficulties of interpretation exist in cuprates, 
where the coherence length is short and the cores may 
nucleate competing order, but in most Fe-based super- 
conductors (FeSC) these complications are less severe or 
absent over a wide range of experimentally tunable pa- 
rameters. 

On the other hand, a complex aspect of these latter 
systems arises due to their multiband nature. The di- 
rectional dependence Vi?(k) also affects the decay length 
of the core states, especially when combined with differ- 
ent gap amplitudes on different Fermi surface sheets. In 
FeSC, the Fermi surface typically consists of two or three 
hole pockets and two electron pockets, as represented in 
the Brillouin zone corresponding to 1-Fe unit cell. The 
size and shape of these pockets varies considerably from 
family to family. A natural question is whether it is the 
normal state band structure and the Fermi surface, or 
the order parameter shape that determine the salient fea- 
tures of the vortex core states as seen in experiment, and 
whether one can draw reliable conclusions about the di- 
rections of the gap nodes or minima based on the real 
space structure of these states. This is the question we 
address in the current Communication. 

The competition between the two effects has been ex- 
plored numerically. For example, the sixfold pattern ob- 
served in 2H-NbSe2 core states [16 can be explained ei- 
ther assuming a weak gap anisotropy or using the angle- 
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FIG. 1. (a) Fermi surface of stoichiometric LiFe As at /c^ = 
in the unfolded 1-Fe "effective" Brihouin zone from DFT. The 
Fermi velocities for different sheets are indicated by the arrows 
pointing to the higher ^(k). We label two inner hole pockets 
ai, a2, one outer hole pocket 7 and two electron pockets 
^2- (b) The Fermi velocity direction 0^^ vs the momentum 
k azimuthal angle for the LiFeAs 7 pocket and the circular 
Fermi surface (shown as insets). 



dependent density of states around the Fermi surface [T7 . 
In pnictides it was argued both that the vortex core states 
are controlled by the order parameter shape [18] and that 
the location of the peak in the DOS is determined by 
the proximity to the band edge in the electron or hole 
bands [19]. To gain qualitative insight into this issue 
we consider a simple model with both the order param- 
eter and band anisotropy characteristic of the Fe-based 
superconductors, and find that in the absence of strong 
nodes the Fermi velocity anisotropy can dominate the 
real-space shape of the vortex core states. 

We focus on the LiFeAs system, which is ideal for STS 
measurements due to its nonpolar surfaces. According 
to density functional theory (DFT) calculations [20 , the 
Fermi surface of this material has three hole pockets and 
two electron pockets, see Fig. [ij The outer hole pocket 
is large and quite square, according to both DFT results 
and ARPES [21 and dHvA [22 measurements. Both 7 
and a2 pocket have small Fermi velocities and therefore 
large normal state DOS. ARPES has identified super- 
conducting leading edge gaps of order 1.5-2 meV for the 
hole pockets, and 3 meV for the electron pockets [21 . 
The London penetration depth data [23 and specific 
heat measurements [24 ruled out the existence of gap 
nodes and were fit to models with two isotropic gaps with 
(Ai,A2) (3meV,1.5meV) and (2 meV, 0.5 meV), re- 
spectively. This suggests moderate gap anisotropy, which 
is not easily detected by the bulk measurements, but can 
substantially affect the real space structure of the core 
states. 

For circular Fermi surfaces the low-energy core bound 
states extend furthest in the direction of the smallest 
gap, but for realistic bands the Fermi velocity anisotropy 
plays a significant role. Since the cross-sections of the 
Pi and P2 electron pockets rotate by a full 180° along 



the kz direction, and since these gaps are larger, it is 
unlikely that these sheets contribute substantially to the 
spatial anisotropy. We therefore focus on the possible 
anisotropy of the gap on the hole Fermi surfaces. The 
most likely candidate for the anisotropic gap that domi- 
nates the low-energy vortex bound states is the 7 pocket. 
The orbital content of this pocket is exclusively dxy^ and 
it couples only weakly to the primarily dxz and dyz elec- 
tron pockets which provide the main pairing weight in 
the conventional spin fluctuation approach 0. It is also 
nearly square, with weakly dispersive parallel surfaces 
oriented along the [110] direction in the 1-Fe zone, and 
with significant variations of the Fermi velocity between 
[100] and [110] directions. Hence we first neglect other 
FS sheets, and contrast the results obtained for the 7 
sheet alone with those for a single circular FS. 

Model. We follow the approach of Ref. [25] that relied 
on the quasiclassical method for superconductivity [26]- 
[28j . used previously to study vortex cores [29^. The 
energy-integrated normal and anomalous Green's func- 
tions ^(r, ^, iuon) and /(r, ^, iuOn) obey the coupled Eilen- 
berger equations 

2 (^ujn + • A(r)^ + ifwF • V /(r, 6>, iujn) 
= 2ig{r,0,iuJn)A{r,0), (la) 



2 (^iujn + • A(r)^ - ihvF • V /(r, 6>, iun) 



2i^(r,6>,icJn)A*(r,6>), 



(lb) 



together with the normalization condition g'^ + // = 1. 
Here A(r) is the vector potential, v^? is the Fermi veloc- 
ity at the location at the Fermi surface labeled by 0, and 
ujn = (2n + l)7rkBT are fermionic Matsubara frequen- 
cies. The Fermi velocity, Vi?(^), is along the 2D unit 
vector k for the circular Fermi surface, and is computed 
for the 7-band in LiFeAs using the Quantum EXPRESSO 
package [30 , as in Ref. 31, In the low field regime, we 
consider the problem of an isolated vortex and assume a 
separable momentum and coordinate dependence of the 
order parameter A(p, k) = Ao$(^) tanh (p/7^^^0), where 
Aq is the bulk gap value in the absence of the field and 
^{0) describes the gap shape on the Fermi surface, = 
l,^d = V2cos26>, and ^s,ani = (1-r cos46>)/vT + 7V2 
with r = 0.3, for the isotropic s-wave, nodal d-wave, and 
extended s-wave gaps respectively. The coherence length 

is ^0 = ^VF,rms/Ao whcrC VF,Tms = \J p (^)\^) F S , 

and the brackets denote the normalized average over the 
Fermi surface, 
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and p{9) is the angle-dependent 



density of states. The factor 77^ accounts for the shrinking 
of core size at low temperature (Kramer- Pesch effect [32l 
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[33]), and we set r]r = 0.1 corresponding to T ~ O.lTc. 
In a fully self-consistent calculation the gap anisotropy 
in momentum space will induce weak core anisotropy in 
real space [34], which we ignore here since the effect is 
small even for nodal systems [34 . 

We solve Eq. ([T]) using the Riccati parameteriza- 
tion [35i and integrating along classical trajectories, 
t{x) = To -\- xvp to obtain the functions g and / at Mat- 
subara frequencies. The LDOS is found after analytic 
continuation from the retarded propagators, N{r^uj) = 
No{Re ^^(ki?, r, + iS))FS' At each point r = (p, (/)) the 
LDOS is obtained by summation over the quasiclassical 
trajectories passing through r. Each trajectory follows 
the direction of the Fermi velocity at a given point on the 
FS, Vi?(k), and samples the gap A(r(x),k). Trajectories 
sampling regions of small order parameter contribute to 
the low energy LDOS. This occurs if the trajectory either 
passes in the vicinity of the core where the order parame- 
ter is suppressed in real space, A(p) <C Aq (small impact 
parameter, dominant for isotropic gaps), or is along the 
direction where the gap has a node or a deep minimum 
in momentum space, A(k) <C Aq (dominant for nodal 
superconductivity) . 

The influence of the FS shape is then clear: the number 
of trajectories with a given impact parameter depends on 
the band structure. Denote the angle between and kx 
axis as 0^^,. For a circular FS = quasiclas- 
sical trajectories in different directions 0^^ are equally 
weighted in FS averaging. In contrast, for anisotropic 
cases, such as the square 7-sheet in LiFeAs, large parts 
of the FS have the along the diagonals (see Fig. [l]3), 
and therefore the average over the trajectories is heavily 
weighted towards that direction as well. 

For an isotropic gap A(k) = const, the largest con- 
tribution to the low energy LDOS at r = (p, 0) comes 
from the trajectories passing through the core, 6>vf = 
or (/) + TT. For a cylindrical FS parameterized by angle 
this corresponds to two points since 0^^, = 0. On an 
anisotropic FS, such as the 7 pocket in LiFeAs, many 
different momentum angles correspond to 6>v^ ^ if , 
and quasiparticles from a large portion of the FS travel 
along these directions. For real space direction (j) = f , all 
these trajectories sample the core region and contribute 
to the low energy LDOS. For ^ away from these directions 
these trajectories have a nonzero impact parameter and 
therefore small weight at low energies. For the extended 
5- wave gap model with r > in the form factor ^s,ani^ 
this implies that the regions of large gap will be empha- 
sized due to preferential directions of Vi?, and therefore 
the FS effects compete with the gap shape in determin- 
ing the spatial profile of the vortex core states. Simply 
assuming that direction of the smallest gap in k space 
yields the orientation of the tails of the bound state wave 
function need not be correct, and may be wrong with a 
strongly anisotropic Fermi surface [? ]. 

Results. Fig. [2] shows the zero energy density of states 




FIG. 2. Normalized ZDOS in a 2.5^o x 2.5^o region around 
the center of the single vortex for different gap models with 
a circular Fermi surface (a - c) and LiFeAs 7 pocket (d - 
f): (a, d) an isotropic s-wave gap Ao; (b, e) a nodal d- 
wave gap Aoa/2cos2^; (c, f) extended s-wave gap Ao(l — 
r cos 4^)/ yl + r^/2, r = 0.3. The gap bulk value is taken to 
be Ao = 1.76Tc. The inset on each panel represents a car- 
toon of the corresponding gap along the Fermi surface. White 
contour lines shown correspond to 0.025A^o. 



(ZDOS) of a circular Fermi surface (a-c) and LiFeAs 7 
pocket (d-f). Comparing panel (a, d) for the isotropic 
gap, we see that the rotation symmetry of ZDOS in (a) 
is broken due to the anisotropy of 7 pocket and Fermi 
velocity; at the same time the ZDOS still preserves the 
crystal four-fold symmetry. In the d-wave case (b) for a 
circular Fermi surface, we recover well-known results for 
the ZDOS, including the double tails along the nodal di- 
rections forced by the vanishing of the bound state wave- 
function exactly along the 45° directions in the quasiclas- 
sical theory [29]. While this feature remains, it becomes 
essentially invisible in the case of the square Fermi surface 
shown in panel (e), as the Fermi surfaces concentrates the 
quasiparticle trajectories even more in the nodal direc- 
tions. Our primary results are now contained in panels 
(c) and (f). The extended-s state ^s,ani has been chosen 
deliberately to have gap minima along the 0° directions 
(along the Fe-Fe bond in the FeSC case). This is clearly 
visible in the case of an isotropic pocket (c), as the tails. 
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while not as well-defined as in the true nodal case, ex- 
tend clearly along these directions in real space. These 
directions rotate by 45°, however, when the same gap 
exists on the square LiFeAs 7 pocket, as in (f). In fact, 
the ZDOS in panel (f) strongly resembles the structure 
observed by Hanaguri et al. in recent STS measurements 
on LiFeAs |36]. 

The results in Fig. [2] strongly challenge the common 
interpretation of STS images of vortices, which assign 
gap minima to the directions of the extended intensity 
in real space. This is probably reasonable in the case of 
true nodes, as indicated by the d-wave examples shown, 
but fails if these minima are not sufficiently deep due to 
the competition with the Fermi surface effects. Now that 
the basic structure of this competition in the case of the 
the ZDOS has been understood, it is interesting to ask 
what may happen in the case of finite energies uj ^ 0. 
The quasiclassical theory incorporates the FS properties 
solely via v^, and thus does not account for the possible 
changes in the shape of the constant energy surfaces for 
STS biases away from zero. Provided the band shape 
varies very slowly on the scale of Tc, this neglect should 
not significantly affect the shape of vortex bound states 
at nonzero energy, however. On the other hand, even 
within the current model, a more important effect may 
be included. In our analysis of LiFeAs, we have until 
now neglected all Fermi surface pockets except the outer 
(7) hole pocket, due to its square shape and because it 
seems likely to have the smallest gap. When the bias 
is increased, higher energy quasiparticle states, including 
those associated with larger gaps, will be probed. Within 
spin fluctuation theory [5 , both the high density of states 
a2 pocket, and the electron pockets, tend to have gap 
minima along the 0° directions. Thus as higher energies 
are probed, it is possible that rotations of the bound state 
shape may take place as the balance between gap struc- 
ture and Fermi surface anisotropy is altered. Unfortu- 
nately even qualitative statements depend on the details 
of the sizes of gaps and gap anisotropics on each sheet, as 
well as on the various Fermi velocities for each band. The 
LiFeAs system is quite clean, however, and if the current 
controversy between ARPES [21 and dHvA [22^ regard- 
ing the Fermi surface can be resolved, spectroscopies of 
bound states on this system should provide enough in- 
formation to determine fairly detailed structure of the 
gap. 

Conclusions. We have used quasiclassical methods to 
calculate the vortex bound states within a single vortex 
approximation, and highlighted the competition between 
gap and Fermi surface anisotropy in the determination 
of the shape of STS images of vortex bound states. If 
the Fermi surface anisotropy is large enough, we have 
shown that the tails of vortex bound states at low energy 



need not correspond to the smallest gaps in the system, 
if those gaps are not true nodes. The ZDOS shape mea- 
sured by STS in experiments on the LiFeAs system with 
very clean surfaces is well reproduced by numerical cal- 
culation. Within our model, we attribute the tail-like 
spectrum to the effect of the non-uniformly distribution 
of Fermi velocity direction on the Fermi surface of the 
LiFeAs 7 hole pocket. Further measurements of the en- 
ergy dependence of bound state shape may further help 
identify the gap anisotropy. 
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